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The region surrounding the center of the Milky Way is both astrophysically rich and complex, 
and is predicted to contain very high densities of dark matter. Utilizing three years of data from the 
Fermi Gamma Ray Space Telescope (and the recently available Pass 7 ultraclean event class), we 
study the morphology and spectrum of the gamma ray emission from this region and find evidence of 
a spatially extended component which peaks at energies between 300 MeV and 10 GeV. We compare 
our results to those reported by other groups and find good agreement. The extended emission could 
potentially originate from either the annihilations of dark matter particles in the inner galaxy, or 
from the collisions of high energy protons that are accelerated by the Milky Way's supermassive 
black hole with gas. If interpreted as dark matter annihilation products, the emission spectrum 
favors dark matter particles with a mass in the range of 7-12 GeV (if annihilating dominantly to 
leptons) or 25-45 GeV (if annihilating dominantly to hadronic final states). The intensity of the 
emission corresponds to a dark matter annihilation cross section consistent with that required to 
generate the observed cosmological abundance in the early universe (cr« ~ 3 X 10"^^ cmVs). We 
also present conservative limits on the dark matter annihilation cross section which are at least as 
stringent as those derived from other observations. 
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I. INTRODUCTION 

Since its launch in June of 2008, the Fermi Gamma 
Ray Space Telescope (FGST) has been producing the 
most detailed and highest resolution observations to date 
of the gamma ray sky between 50 MeV and 100 GeV. 
Among the objectives of this experiment are to increase 
our understanding of how astrophysical objects such as 
black holes and pulsars accelerate cosmic rays, and to 
help identify the substance or substances that compose 
the dark matter of our universe. For each of these areas 
of inquiry, the region surrounding the center of the Milky 
Way represents a particularly interesting and promising 
target of study. On the one hand, the Galactic Center 
is an extraordinarily rich and complex region, containing 
our galaxy's supermassive black hole, as well as super- 
nova remnants, massive X-ray binary systems, massive O 
and B type stars, and two young and massive star clus- 
ters (Arches and Quintuplet) [IHl]. On the other hand, 
the Galactic Center is predicted to contain very high den- 
sities of dark matter, which in many models leads to a 
very high rate of dark matter annihilation, and a cor- 
respondingly high luminosity of gamma rays. No other 
astrophysical source or region is expected to be as bright 
in dark matter annihilation products as the Galactic Cen- 
ter. 

In this article, we follow previous work [5 and perform 
a detailed study of the spectral and morphological fea- 
tures of the gamma rays from the Galactic Center region, 
with the intention of identifying or constraining the ori- 
gins of these particles. In particular, we produce gamma 
ray maps which reveal the presence of both a bright, ap- 
proximately point-like, gamma ray source at the Galactic 
Center, along with a more spatially extended emission 



component. The spectrum of this extended source peaks 
strongly between several hundred MeV and ^10 GeV. 
We find good agreement between our results and those 
reported by other groups (SHZI- 

In discussing the possible origins of this extended emis- 
sion, we find again that the observed spectrum and mor- 
phology are consistent with that predicted from annihi- 
lating dark matter particles with a mass of 7-12 GeV 
annihilating dominantly to leptons [5] or a mass of 25-45 
GeV annihilating dominantly to hadronic final states [H] . 
In either case, the normalization of the gamma ray flux 
requires an annihilation cross section that is consistent, 
within astrophysical uncertainties, with the value pre- 
dicted for a simple thermal relic {av ^ 3 x 10^^^ cm^/s). 
We also discuss the possibility that the extended gamma 
ray emission is produced through the collisions of ener- 
getic protons which are accelerated by the supermassive 
black hole with gas [6 . While we consider this to be 
the leading astrophysical explanation for the gamma ray 
emission observed by the FGST, it is somewhat difficult 
to assess this hypothesis given how little is known or can 
be reliably predicted about the spectrum or flux of pro- 
tons accelerated by the central black hole, and how little 
is known about the history of this object (such as peri- 
ods of flaring and relative inactivity) and the properties 
of the surrounding interstellar medium. 

The remainder of this article is structured as follows. 
In Sec. [nj we describe our analysis of the Fermi data and 
present gamma ray maps of the Inner Galaxy and the cor- 
responding spectrum of this emission. In Sec. |III| we fur- 
ther describe the properties of this emission and compare 
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our results to those found by other groups. In Sec 
we discuss several possible origins of this emission, in- 
cluding energetic protons from the central supermassive 



2 



black hole, dark matter annihilations, and a population 
of gamma ray pulsars. In Sec. |Vj we derive constraints 
on the dark matter annihilation cross section which are 
at least as stringent as those based on other observa- 
tions, such as those of dwarf spheroidals, galaxy clusters, 
the cosmological diffuse background, and nearby subha- 
los. In Sec. |VI[ we discuss our results within the larger 
context of dark matter searches and summarize our con- 
clusions. 



II. ANALYSIS PROCEDURE 

We begin our analysis by generating contour maps of 
the region surrounding the Galactic Center which de- 
scribe the distribution of gamma rays observed by the 
Fermi-LAT (Large Area Telescope) over the three years 
between August 4, 2008 and August 3, 2011. These 
maps were derived using only front-converting events 
(which have a superior point-spread function compared 
to back-converting events) from the Pass 7 ultraclean 
class. As recommended by the FGST collaboration, we 
include only events with zenith angles smaller than 100 
degrees, and do not include events recorded while the 
Fermi satellite was transitioning through the South At- 
lantic Anomaly or while the instrument was not in survey 
mode. Each of the maps has been smoothed out at a scale 
of 0.5 degrees (the contour maps thus represent the flux 
observed within a 0.5 degree radius of a given direction in 
the sky). These raw maps are shown in the left frames of 
Fig. [T] for five different energy ranges between 100 MeV 
and 100 GeV. 

In each map, ten contours are shown, distributed 
linearly between 2.64 x 10~^ and 2.64 x 10^''cm^^s~^ 
sq deg-i (100-300 MeV), 2.45 x 10"^ and 2.45 x 
lO-'^cm-^s-i sqdeg-i (300-1000 MeV), 1.07 x 10"^ and 
1.07 X 10-^cm-2s-i sq deg-^ (1-3 GeV), 2.66 x 10"^ 
and 2.66 x lO-^cm-^s-^ sq deg-^ (3-10 GeV), and 
3.77 X 10-10 and 3.77 x lO-'^cm-^g-i sq deg-^ (10-100 
GeV). Note that the 2.64 x lO^^cm^^ s^^ sq deg"^ con- 
tour appears out of the field in the upper-left and upper- 
middle frames. 

The blue points shown in the maps represent the lo- 
cations of sources contained in the Fermi Second Source 
Catalog (2FGL) [3], and the size of each point is propor- 
tional to the reported intensity of the source in the energy 
range shown. To account for these sources, we have gen- 
erated a template map of their emission (assuming the 
central values for their intensity and locations as reported 
in the 2FGL), and taking into account the point-spread 
function of the Fermi-LAT (as determined by the Fermi 
Tool gtpsf). In the center frames of Fig. l] we show the 
maps as they appear after subtracting this source tem- 
plate. Note that we have not removed the central bright 
source, as its emission is difficult to disentangle from dark 
matter annihilation products originating from the inner 
region of a cusped halo profile. We will return to this 
issue later in the article. 



After subtracting these known sources, the dominant 
remaining component is the diffuse emission associated 
with the disk of our galaxy. This emission is dominated 
by cosmic ray processes taking place throughout the disk 
of the Milky Way, which one must look through in order 
to observe the Galactic Center. By studying the mor- 
phology of this emission over the regions of 5° < |Z| < 10°, 
we find only a modest degree of variation with galactic 
longitude. In Fig. [2] we show as solid lines the observed 
gamma ray flux as a function of galactic latitude, aver- 
aged over the range of 7° < |/| < 10° (in order to avoid 
any contamination with emission from the inner most 
degrees, we do not here make use of the data within 7°). 

The gamma ray emission from the disk of our galaxy 
is dominated by the decays of neutral pions produced in 
cosmic ray interactions with gas, although inverse comp- 
ton and bremsstrahlung components also contribute. To 
model the morphology of the pion decay component, we 
adopt the following distribution of gas [lOl [11] : 



Pgas 
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R < 7kpc, 
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Pgas oc e-l^l/^=^(«)e-«/«-, i?>7kpc, 



where z and R describe the location relative to the Galac- 
tic Center in cylindrical coordinates. We set Rgc = 3.15 
kpc (as fit to the data shown in Fig. 4 of Ref. [TUj) and 
z^^{R) = 0.1 + 0.00208 X (i?/kpc)2 kpc (as fit to Fig. 4 of 
Ref. [H]), in good agreement with observations of 21-cm 
surveys, which trace the density of neutral hydrogen. To 
estimate the flux of pion decay gamma rays, we integrate 
this distribution over the line-of-sight (and again smooth 
over a radius of 0.5 degrees). After accounting for the 
Fermi-LAT point-spread function, we find that this gas 
distribution leads to the morphology described by the 
dashed lines in Fig. [2] This is in good agreement with 
the observed morphology of the diffuse emission. We also 
note that the spectral shape implied by the relative fluxes 
in these five energy bins is consistent with that predicted 
for a combination of pion decay and inverse compton 
scattering processes, as previously found in Ref. [S]. By 
subtracting this disk template from the gamma ray maps, 
we are able to remove the overwhelming majority of the 
diffuse astrophysical background from our maps. We em- 
phasize that in performing this subtraction, we are not 
extrapolating any physical features of the inner galaxy, 
but are merely extrapolating the line-of-sight gas densi- 
ties along the disk from directions slightly away from the 
Galactic Center to those more aligned with the Galactice 
Center. 

In the right frames of Fig. [T] we show the resulting 
maps after subtracting both the known sources template 
(again, not including the bright central source) and the 
line-of-sight gas template. In each energy range, the ma- 
jority of the background has been accurately removed 
by this simple subtraction. While this subtraction pro- 
cedure does not perfectly remove all likely astrophysical 
backgrounds, the residuals outside of the inner ~2° are 
very modest, typically on order of 10% or less of the 
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FIG. 1: Contour maps of the gamma ray flux from the region surrounding the Galactic Center, as observed by the Fermi 
Gamma Ray Space Telescope. The left frames show the raw maps, while the center and right frames show the maps after 
subtracting known sources (not including the central source), and known sources plus emission from cosmic ray interactions 
with gas in the Galactic Disk, respectively. All maps have been smoothed over a scale of 0.5 degrees. See text for more details. 



residual flux in the innermost region of the Galaxy. We 
include the observed spatial variations of the residuals as 
a systematic error, which we propagate throughout this 
study. 

The residuals in this innermost region include 
a roughly spherically symmetric component centered 
around the Galactic Center, along with a sub-dominant 
component that is somewhat extended along the disk. 
Due to its similar angular extent, we consider it likely 
that this component is associated with emission from 
proton-proton collisions taking place in the Galactic 



Ridge, as observed at higher energies by HESS [12]. The 
remaining spherically symmetric component could plau- 
sibly originate from dark matter annihilations, processes 
associated with the Milky Way's supermassive black hole, 
gamma ray pulsars, or a combination of these and other 
sources. We will return to these issues in Sees. IIIII and 

In Fig. [Sj we show the spectrum of the emission from 
the inner 5 degrees surrounding the Galactic Center, after 
removing the known sources and disk emission templates. 
The spectrum is clearly brightest between 300 MeV and 
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FIG. 2: The observed gamma ray flux (after subtracting known sources) averaged over 7° < \l\ < 10° as a function of galactic 
latitude (solid), compared to that predicted from the line-of-sight gas density (dashes). See text for details. 



10 GeV, and drops by nearly an order of magnitude above 
^10 GeV. Note that the spectral shape of this residual is 
very similar to that (preliminarily) reported in conference 
presentations by the Fermi Collaboration [13]. In the 
following sections, we will explore the characteristics of 
this residual emission and discuss its possible origins. 

III. PROPERTIES OF THE INNER EMISSION 

In an effort to constrain the origin (or origins) of the 
gamma rays from the inner region of our galaxy, we dis- 
cuss in this section the spectral and morphological char- 
acteristics of the observed emission [14]. We begin by 
noting that the morphology of the observed residual is 
not consistent with that of a single point source. In par- 
ticular, we find that above 300 MeV, less than half of 



the residual emission shown in the right frames of Fig. [T] 
can be accounted for by a single, centrally-located point 
source. This conclusion is also supported by the inde- 
pendent analyses by Boyarsky et al. [7], Chernyakova 
et al. fS^, and Hooper and Goodenough which each 
found a spectrum of point-like emission from the Galactic 
Center which is considerably less intense than the total 
residual emission shown in Fig. [3] In Fig. [4] we show the 
spectra of point-like emission from the Galactic Center, 
as reported in each of these three prior studies. We note 
that the intensity and spectral features of the Galactic 
Center point source found by these three groups are very 
similar, despite the very different analysis techniques em- 
ployed. 

In Fig. [5] we compare this spectrum of point-like emis- 
sion (as reported by Boyarsky et al. [7]) to the spectrum 



5 




6 

o 

% 10" 
\ 



(GeV) 

FIG. 3: The spectrum of the residual emission from the inner 
5 degrees surrounding the Galactic Center, after subtracting 
the known sources and line-of-sight gas templates. 
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FIG. 4: Fits for the spectrum of the central emission, assum- 
ing a point-like source morphology, from the previous work 
of three different groups [SHZj. Despite the different analysis 
approaches taken, these fits are all in reasonable agreement. 
The dashed line is the broken power-law fit to this spectrum 
as presented in Ref. [6]. 
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FIG. 5: A comparison of the total residual emission found in 
this study (black) with the spectra of point-like emission (red) 
and extended emission (blue) (as in the case of annihilating 



dark matter with pDM oc r 



as presented in Ref. [7]. This 



comparison supports our finding that this residual emission 
below ~ 300 MeV is consistent with a point-like source ori- 
gin, while much of the emission at higher energies is indeed 
spatially extended. 



tially extended component in their model (with a mor- 
phology corresponding to that predicted for annihilating 
dark matter with a distribution given by poM oc r^^ '^*), 
they found that the fit improved considerably (reducing 
the log-likelihood by 25 with the addition of only one 
new parameter) [7]. The spectrum of this spatially ex- 
tended component is also shown in Fig. [5] The spectrum 
of the residual emission found in our analysis is in very 
good agreement with the sum of point-like and extended 
components as reported by Boyarsky et al. From these 
comparisons, we conclude that in addition to the presence 
of point-like emission from the Galactic Center, a com- 
ponent of extended emission is also prominently present 
at energies greater than ^300 MeV.^ 



of residual emission found in our analysis.^ Between 100- 
300 MeV, there is good agreement, indicating that most 
or even all of the residual gamma rays in this energy range 
could originate from a single point source. At higher 
energies, however, the residual emission consistently ex- 
ceeds the flux attributable to point-like emission; by a 
factor of ~2-3 between 0.3 and 3 GeV, and by a factor of 
^5 above 3 GeV. When Boyarsky et al. included a spa- 



IV. POSSIBLE ORIGINS OF THE OBSERVED 
EMISSION 

A number of proposals have been put forth to explain 
the bright gamma ray emission observed from the Galac- 
tic Center by the FGST. These possibilities include the 
central supermassive black hole tS, a population of 



HESS |15l and other ground based telescopes |16| have also ob- 
served point-like emission from the Galactic Center at energies 
above ~200 GeV. This very high energy gamma ray source may 
be associated with the point-like emission observed at lower en- 
ergies, as shown in Fig. [4| 



^ The spectrum of our residual as presented in Figs.[3|and[5]denotes 
the residual within a 5 degree radius around the Galactic Center, 
whereas the spectrum of extended emission reported in Ref. [7| 
is taken from a similar, but not identical, inner 10° X 10° region. 
Given the highly concentrated nature of the morphology being 
considered, however, this difference is negligible. 
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unresolved millisecond pulsars [T7], and dark matter an- 
nihilations [3 [8]. In this section, we discuss the mor- 
phological and spectral characteristics of the gamma ray 
emission expected from each of these potential sources 
and compare this to the emission observed by the FGST. 



A. Cosmic ray Acceleration by the Supermassive 
Black Hole 

To begin, we reiterate that the morphology of the ob- 
served emission is not entirely point-like in nature, but 
instead is somewhat spatially extended. This allows us to 
rule out the possibility that gamma rays directly emitted 
by the Milky Way's central black hole are responsible for 
the observed emission.'^ If, however, the observed gamma 
ray spectrum originates from cosmic rays that have been 
accelerated by the black hole, then a spatially extended 
distribution of gamma rays could result. 

For example, it was previously proposed that the TeV- 
scale gamma rays observed from the Galactic Center 
could originate from the inverse Compton scattering of 
energetic electrons accelerated by the black hole [111 [5D] . 
This scenario, however, predicts considerably less GeV- 
scale emission than is observed by Fermi, and thus can- 
not account for the residual emission discussed here [5]. 
Alternatively, the black hole may accelerate cosmic ray 
protons which then diffuse throughout the surrounding 
interstellar medium, producing pions and thus gamma 
rays through interactions with gas [51 [H]. The spec- 
trum and spatial distribution of the gamma ray emis- 
sion resulting from this process depends not only on the 
spectrum of protons injected from the black hole, but 
also on the diffusion coefficient and distribution of gas 
in the surrounding medium, as well as on the emission 
history of the black hole (occurrences of flares and peri- 
ods of relative dormancy). As none of these inputs are 
currently very well constrained, it is difficult to make re- 
liable predictions for the resulting gamma ray spectrum 
and distribution. That being said, it appears plausible 
that a reasonable astrophysical scenario could potentially 
explain much of the observed emission. 

Perhaps the greatest challenge in accounting for the 
observed emission with energetic protons accelerated by 
the central black hole is the very rapid increase in the flux 
of spatially extended emission observed between approx- 
imately 200 and 700 MeV (see blue error bars in Fig. [5]). 
Even for a mono-energetic spectrum of protons, the re- 
suling spectrum of gamma rays from pion decay does not 
rise rapidly enough to account for this feature. Perhaps 



^ Further supporting this conclusion is the fact that the cen- 
tral emission observed by Fermi shows no variability on month 
timescales [6], as would be expected based on the variability of 
this source in X-ray and infrared wavelengths 1181 . The emis- 
sion from a ^SO parsec source region would have any variability 
suppressed on timescales shorter than Rjc ~ 100 years. 



this could be reconciled, however, if a sizeable fraction 
of the apparently point-like emission in the 100-300 MeV 
bin is in fact somewhat extended and arises from cosmic 
ray interactions. 

Lastly, we also note that a sizable fraction of the high 
energy emission observed by the FGST is likely to be as- 
sociated with the HESS galactic ridge. This ridge emis- 
sion, as measured by HESS, possesses a power- law-like 
spectrum with a spectral index of 2.29±0.07stat ±0.20sys 
over the energy range of approximately 0.2 to 10 TeV. 
Due to the spatial correlation of this emission with the 
locations of molecular clouds in the central 200 parsecs 
of the Milky Way, the origin of the ridge emission is con- 
ventionally taken to be the decays of neutral pions pro- 
duced in the interactions of cosmic ray protons or nuclei 
with the surrounding molecular gas. In order to generate 
a gamma ray spectrum with this spectral index, the re- 
sponsible protons are required to have a spectral index of 
approximately 1.9 ±0.2 [5^. Based on an extrapolation of 
this spectral shape, we estimate that on the order of 30% 
of the 10-100 GeV residual is associated with the ridge. 
At energies below ~10 GeV, however, the ridge emis- 
sion constitutes a much smaller fraction of the observed 
residual, unless the spectrum of cosmic ray protons in the 
region is significantly enhanced below ~50 GeV relative 
to the power-law behavior we have assumed. 

B. Annihilating Dark Matter 

It has long been appreciated that if dark matter par- 
ticles annihilate in pairs (as predicted in most models 
of weakly interacting massive particles), the resulting 
gamma ray signal would be brightest from the direction 
of the Galactic Center [22]. The energy and angular de- 
pendent flux of such gamma rays is given by 

where (crw) is the dark matter annihilation cross section 
multiplied by the relative velocity of the two dark matter 
particles (averaged over the velocity distribution), moM 
is the mass of the dark matter particle, i^) is the angle 
observed relative to the direction of the Galactic Center, 
p(t) is the dark matter density as a function of distance 
to the Galactic Center, and the integral is performed over 
the line-of-sight. dN^/dE^ is the gamma ray spectrum 
generated per annihilation, which depends on the mass 
and dominant annihilation channels of the dark matter 
particle (we use PYTHIA 23J to calculate dN^/dEj for 
various dark matter scenarios in this study). 

Modern numerical simulations of the evolution of cold 
dark matter predict the formation of halos with a nearly 
universal density profile [24] . Within the inner volumes 
of such halos, the density of dark matter varies as puM oc 
r^'' , where r is the distance to the halo's center. The 
frequently used Navarro Frenk and White (NFW) profile, 
for example, features an inner slope of 7 = 1.0 j25j . The 
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results of the Via Lactea II simulation favor a somewhat 
steeper inner slope (7 « 1.2) [26], while the Aquarius 
Project finds a somewhat less steep value which varies 
with r [37]. 

When considering the dark matter distribution in the 
central kiloparsecs of the Milky Way, it is important to 
include the effects of stars and gas, which are not taken 
into account by dark matter-only simulations such as Via 
Lactea II and Aquarius, but which dominate the gravita- 
tional potential of the Inner Galaxy. Generally speaking, 
as a result of dissipating baryons, dark matter density 
profiles are expected to be adiabatically contracted, re- 
sulting in the steepening of their inner profiles [28]. The 
degree to which this effect is manifest depends on the 
fraction of the baryons that dissipate slowly by radiative 
cooling. 

As hydrodynamical simulations which model the pro- 
cess of galaxy formation have improved, efforts to predict 
the effects of baryons have begun to converge. In partic- 
ular, several groups (using different codes) have consis- 
tently found that Milky Way sized halos are adiabatically 
contracted, increasing the density of dark matter in their 
inner volumes relative to that predicted by dark matter- 
only simulations (see Ref. [S^ and references therein). 
These simulations, which include the effects of gas cool- 
ing, star formation, and stellar feedback, predict a de- 
gree of adiabatic contraction which steepens the inner 
slopes of dark matter density profiles from 7 « 1.0 to 
7 « 1.2 — 1.5 within the inner ^10 kpc of Milky Way-like 
galaxies [IHl [50] . The resolution of such simulations is 
currently limited to scales larger than ^100 parsecs [5T] . 

With this information in mind, we can compare the 
expected spatial distribution of dark matter to the ob- 
served angular distribution of gamma rays from around 
the Galactic Center. Making this comparison, we find 
that the majority of the residual emission observed be- 
tween 300 MeV and 10 GeV can be described by annihi- 
lating dark matter with a distribution given by p{r) cx r'^ , 
with 7 « 1.25 — 1.40.^ In contrast, an NFW-like profile 
with 7 = 1.0 would predict a considerably broader distri- 
bution of gamma rays than is found in our residual maps. 
More quantitatively speaking, for 7 = 1.0 we find that 
for energies of 300-1000 MeV, 1-3 GeV, and 3-10 GeV, 
respectively, no more than 22%, 18% and 27% of the flux 
found in the innermost half degree around the Galactic 
Center can arise from dark matter annihilations without 
also exceeding the flux observed at distances beyond one 
degree. In contrast, if we select an inner slope of 7 = 1.3, 
we find that up to 72%, 74% and 100% of the innermost 
emission could originate from dark matter annihilations. 
The remainder of the residual could easily originate from 
the central point source with the spectrum presented in 



* Below 300 MeV, the observed emission is dominated by point- 
like emission, and the flux of the emission drops off significantly 
above 10 GeV, leading us to focus on this energy range. 
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FIG. 6: The range of dark matter masses and annihilation 
cross sections for which dark matter annihilations can account 
for the majority of the observed residual emission between 
300 MeV and 10 GeV, for three choices of annihilation chan- 
nels ("leptons" denotes equal fractions to e'^e~, fi^ n~ and 
T^T~). Also shown for comparison is the annihilation cross 
section predicted for a simple thermal relic {av = 3 x lO"'^^ 
cm^/s). Note that there is a factor of a few uncertainty in the 
annihilation cross section, corresponding to the overall dark 
matter density and distribution. See text for details. 



Fig. [4] 

If a sizable fraction of the residual emission does origi- 
nate from annihilating dark matter, then we can use the 
spectrum of this emission to inform us as to the mass and 
dominant annihilation channels of the dark matter par- 
ticles. In particular, the rapid decrease in the flux above 
^10 GeV suggests that the spectrum is being dominated 
by ~30 GeV dark matter particles annihilating to quarks, 
or by ^10 GeV particles annihilating to leptons (among 
annihilations to leptons, those to taus produce far more 
gamma rays than those to either muons or electrons). 
In Fig. [6] we show the range of dark matter masses and 
annihilation cross sections for which dark matter anni- 
hilations can account for the majority of the observed 
residual emission (without exceeding the observed resid- 
ual) in each of the three energy bins between 300 MeV 
and 10 GeV, for three choices of the annihilation chan- 
nels. Interestingly, we note that the normalization of the 
signal requires us to consider annihilation cross sections 
that are within a factor of a few of the value predicted 
for a simple thermal relic {av = 3 x 10^^^ cm'^/s). The 
precise value of the required annihilation cross section 
depends on the quantity of dark matter present, and is 
thus subject to the related uncertainties. In Fig. [6] and 
throughout the remainder of this paper, we have normal- 
ized the dark matter distribution such that the total mass 
of dark matter within the solar circle is 3.76 x 10^^ GeV, 
which is the value corresponding to the case of 7 = 1.0 
and a local density of 0.4 GeV/cm^. This value is sup- 
ported by a combination of microlensing and dynami- 
cal constraints, although uncertainties exists [35]. With 
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FIG. 7: Examples illustrating how dark matter annihilations and astrophysical sources could combine to make up the observed 
residual emission surrounding the Galactic Center. In the upper left frame, we show results for a 10 GeV dark matter particle 
with an annihilation cross section of ctw = 7 x 10~^^ cm'^/a and which annihilates only to leptons (e^e~, fJ-'^ fJ.~ and r^r~, 
1/3 of the time to each). In the upper right frame, we show the same case, but with 10% of the annihilations proceeding to 
bb. In the lower frame, we show results for a 30 GeV dark matter particle annihilating to bb with an annihilation cross section 
of cru = 6 X 10~^^ cm^/s. In each case, the annihilation rate is normalized to a halo profile with 7 = 1.3. The point source 
spectrum is taken as the broken power-law shown in Fig. |4] and the Galactic Ridge emission has been extrapolated from the 
higher energy spectrum reported by HESS 12 , assuming a pion decay origin and a power-law proton spectrum. See text for 
details. 



these uncertainties in mind, one should consider all an- 
nihilation cross sections shown in Fig. [6] and elsewhere in 
this paper to be accurate only to within a factor of a few. 

Of course, it is also expected that astrophysical sources 
will contribute to the Galactic Center's gamma ray spec- 
trum between 300 MeV and 10 GeV. In Fig. [7] we show 
three examples in which emission from a central point 
source (as shown in Fig.|4]), along with emission from the 
Galactic Ridge (as extrapolated from the higher energy 
HESS emission, assuming a spectral shape that results 
from a power-law spectrum of protons) combine with a 
contribution from dark matter to generate the observed 
residual emission. Note that the lowest energy emission 
is largely generated by the central point source (as sug- 
gested by the observed morphology) while the highest 
energy bin is dominated by emission from the Galactic 
Ridge. Only the 300 MeV-10 GeV range is dominated by 
dark matter annihilation products. 



C. Millisecond Pulsars 

A population of gamma ray point sources surround- 
ing the Galactic Center could also potentially contribute 
to the observed residual emission. Millisecond pulsars, 
which are observed to produce spectra that fall off rapidly 
above a few GeV, represent such a possibility [51 [T7]. 

Observations of resolved millisecond pulsars by FGST 
have found an average spectrum well described by 
dN^/dE^ cx E-^-^ exp{-E^/2.8 GeV) [33 . Similarly, the 
46 gamma ray pulsars (millisecond and otherwise) in the 
FGST's first pulsar catalog have a distribution of spec- 
tral indices which peaks strongly at F =1.38, with 44 
out of 46 of the observed pulsars possessing (central val- 
ues of their) spectral indices greater than 1.0 [34] (see 
Fig. [s]). In contrast, to produce a sizable fraction of the 
spatially extended residual emission between 300 MeV 
and 10 GeV without exceeding the emission observed be- 
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FIG. 8: A histogram showing the distribution of spectral in- 
dices, r, of pulsars in the Fermi Pulsar Catalog. 
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FIG. 9: The spectral indices (with statistical and systematic 
error bars) of the eight globular clusters observed by the Fermi 
Gamma Ray Space Telescope [35] . 



low 300 MeV, the average pulsar in the Galactic Center 
population would be required to possess a spectral index 
harder than F 1.0. And although we agree with the au- 
thor of Ref. [17 that a small number of pulsars (including 
J1958-F2846, J2032+4127 and J2043+2740) have been 
observed with such hard spectral indices, we do not be- 
lieve that the existing data supports the conclusion that a 
large population of pulsars (as would be required to gen- 
erate the observed emission) would produce an average 
gamma ray flux with a spectral shape able to account for 
the observed emission from the Galactic Center.^ That 
being said, if the population of pulsars present in the 
central stellar cluster were to differ significantly from the 
sample represented by the Fermi pulsar catalog, a differ- 
ent conclusion could potentially be reached. 

An opportunity to measure the emission from large 
populations of gamma ray pulsars exists in the form of 
globular clusters, whose gamma ray emission is generally 
attributed to pulsars contained within their volumes. Un- 
fortunately, the gamma ray spectra of these objects have 
not been well measured. In particular, the eight globular 
clusters with spectra reported by Fermi have an average 
spectral index very close to that of pulsars (F w 1.38), 
but with very large individual error bars which extend 
from roughly to 2.5 (these values, including la statis- 
tical and systematic errors are shown in Fig. |9]). Per- 
haps with more data, we will learn from these systems 
whether the spectral indices of large pulsar populations 
can be hard enough to accommodate the emission ob- 
served from the Galactic Center. 

Lastly, we note that it is somewhat difficult to accom- 



^ The error bars on the spectral indices of these three hardest 
pulsars are also quite large, F = 0.77 it 0.31, 0.68 it 0.46, and 
1.07 it 0.66 [34j. 



modate the very spatially concentrated morphology of 
the observed gamma ray emission with pulsars. As origi- 
nally pointed out in Ref. [5j , to match the observed angu- 
lar distribution of this signal, the number density of pul- 
sars would have to fall off with the distance to the Galac- 
tic Center at least as rapidly as r^^'^. In contrast, within 
the innermost parsec of the Galactic Center, the stellar 
density has been observed to fall off only about half as 
rapidly, r~^'^^ |36j . Furthermore, even modest pulsar 
kicks of ~ 100 km/s would allow a pulsar 10 pc from 
the Galactic Center to escape the region, consequently 
broadening the angular width of the signal. Annihilat- 
ing dark matter, in contrast, produces a fiux of gamma 
rays that scales with its density squared, and thus can 
much more easily account for the high concentration of 
the observed signal. 



V. CONSTRAINTS ON THE DARK MATTER 
ANNIHILATION CROSS SECTION 

In this section, instead of attempting to determine the 
origin of the gamma rays from the Galactic Center re- 
gion, we use the observed spectrum and flux to place 
limits on the dark matter annihilation cross section. In 
doing this, we do not assume anything about the source 
or sources responsible for the observed emission, but in- 
stead only require that dark matter annihilation prod- 
ucts do not exceed the observed emission (after subtract- 
ing the known sources and line-of-sight gas templates, 
as described in Sec. [ll|. Despite using this very simple 
and conservative approach, we derive constraints that are 
competitive with or stronger than those placed by other 
indirect search strategies, including those from observa- 
tions of dwarf spheroidals (37|, galaxy clusters |38|, the 
cosmological diffuse background '39J, and nearby subha- 
los go]. 
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FIG. 10: Upper limits on the dark matter annihilation cross section for several choices of the final state and for three values 
of the halo profile's inner slope, 7. Also shown for comparison is the annihilation cross section predicted for a simple thermal 
relic (av = 3 x 10~^® cm^/s). Uncertainties in the overall dark matter density have not been included, but based on the errors 
presented in Ref. [32] . we expect that this would only weaken our limits by about 30-50%. See text for more details. 



In Fig. 10 we show the 95% confidence level upper 
limits on the dark matter annihilation cross section for 
several choice of the final state, and for three values of the 
halo profile's inner slope. Based on the results of hydro- 
dynamical simulations [211130], we consider the j — 1.2 
to represent the minimal degree of baryonic contraction, 
where as the 7 = 1.3 and 1.4 cases should be taken as 
more central estimates. 

The constraints we have derived from the Galactic 
Center region are indeed quite stringent. Even in the 
case of only a very modest degree of baryonic contrac- 
tion (7 = 1.2), we find that dark matter particles with 
the canonical annihilation cross section of av ~ 3x 10^^^ 
cm'^/s and which proceed to hadronic final states are pre- 
dicted to exceed the observed gamma ray fiux from the 
Galactic Genter unless they are more massive than ap- 
proximately 300 GeV. In comparison, the Fermi collab- 



oration's combined analysis of 10 dwarf spheroidals only 
excludes such dark matter particles with masses below 
approximately 30 GeV [37]. ^ 



In Fig. [TT] we also show the constraints which result 
if the effects of baryons on the dark matter distribution 
are neglected (using a uncontracted NFW profile). Even 
in this case, we find limits which are approximately as 
stringent as those derived from dwarf galaxies. 



^ Unlike the central regions of Milky Way-like halos, the dark 
matter density profiles of dwarf spheroidal galaxies are not gen- 
erally expected to be contracted by baryons 1411 . The uncon- 
tracted NFW profile adopted in the Fermi Collaboration dwarf 
spheroidal analysis is thus appropriate. 
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FIG. 11: Upper limits on the dark matter annihilation cross section for several choice of the final state, neglecting the effects 
of baryons (using an uncontracted NFW halo profile). Also shown for comparison is the annihilation cross section predicted 
for a simple thermal relic {av — 3 x lO"'^^ cm^/s). Uncertainties in the overall dark matter density have not been included, 
but based on the errors presented in Ref. |32], we expect that this would only weaken our limits by about 30-50%. See text for 
more details. 



VI. DISCUSSION AND SUMMARY 

In this article, we have used the first three years of 
data taken by the Fermi Gamma Ray Space Telescope 
(FGST) to study the spectrum and spatial morphology 
of the gamma ray emission from the region surround- 
ing the Galactic Center. In doing so, we have identified 
a spatially extended component of gamma ray emission 
which peaks at energies between approximately 300 MeV 
and 10 GeV. The origin of these gamma rays is currently 
uncertain, although they could potentially be the annihi- 
lation products of dark matter particles, or the products 
of collisions of high energy protons accelerated by the 
Milky Way's supermassive black hole with gas. 

If this extended source of gamma rays is interpreted 
as dark matter annihilation products, the spectrum of 
this emission favors dark matter particles with a mass 
in the range of 7-12 GeV (if annihilating dominantly 
to leptons) or 25-45 GeV (if annihilating dominantly to 
hadronic final states). The former of these mass ranges is 
of particular interest in light of the observations reported 
by the direct detection experiments DAMA/LIBRA WI], 
CoGeNT [43;, and CRESST 44, which each report 
signals consistent with an approximately 10 GeV dark 
matter particle (see also, however, constraints from the 
CDMS [15] and XENON gg] collaborations, and related 
discussions [17]). Further motivating the dark matter in- 
terpretation of the Galactic Center gamma rays is the 
fact that the annihilation cross section required to nor- 
malize the annihilation rate to the observed flux is ap- 
proximately equal to the value required to generate the 
observed cosmological abundance in the early universe 
{av ~ 3 X 10~^^ cm'^/s). In other words, in lieu of res- 
onances, coannihilations, P-wave suppression, or other 
complicating factors, a particle species that will freeze- 
out in the early universe with a density equal to the mea- 
sured dark matter abundance is also predicted to anni- 
hilate today at a rate that is similar to that needed to 
produce the observed gamma rays from the Galactic Cen- 



ter. 

Additionally, we point out that if dark matter par- 
ticles are annihilating in the Inner Galaxy at the rate 
required to produce the observed gamma ray flux, then 
the resulting energetic electrons and positrons will dif- 
fuse outward, potentially producing observable quanti- 
ties of synchrotron emission. In particular, focusing on 
the case of 7-12 GeV dark matter particles annihilating 
dominantly to leptons, the halo profile and cross section 
required to produce the morphology and normalization 
of the observed gamma ray flux is also predicted to lead 
to the production of a diffuse haze of synchrotron emis- 
sion very similar to that observed by WMAP [IS] (see 
Fig. 3 of Ref. [49 for a direct comparison). It also ap- 
pears that the excess radio emission observed at higher 
galactic longitudes by the ARCADE 2 experiment [SD] 
possesses a spectral shape and overall intensity consis- 
tent with originating from dark matter with the same 
mass, cross section, dominant channels, and distribu- 
tion [51] [52]. Lastly, we mention that 7-12 GeV dark 
matter particles with the distribution and annihilation 
cross section favored here would be capable of deposit- 
ing the required energetic electrons into the Milky Way's 
non-thermal radio filaments [53 , providing an explana- 
tion for their peculiar spectral features. 

It is noteworthy that the different explanations pro- 
posed for the observed gamma ray emission from the 
Galactic Center predict different accompanying spectra 
of cosmic ray electrons, potentially providing us with a 
way to discriminate between these different scenarios. Of 
the sources proposed for the observed gamma ray emis- 
sion, only dark matter annihilations are predicted to pro- 
duce comparable fluxes of gamma rays and electrons, 
with spectra that peak at similar energies. Pulsars, in 
contrast, produce gamma ray spectra which peak at ~1- 
3 GeV and electron spectra which peak at several hun- 
dred GeV [54]. Perhaps future observations of the Inner 
Galaxy at radio and microwave frequencies will be able 
to make use of this comparison to shed light on the origin 
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of the gamma ray emission from the center of our galaxy. 

Lastly, we have also presented conservative limits on 
the dark matter's annihilation cross section which are at 
least as stringent as those derived from other observa- 
tions, such as those of dwarf spheroidal galaxies. 
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